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^  Experiments  showing  improved  levels  of  population  inversion  on  soft  x-ray 
transitions  in  laser  heated  plasmas  were  performed  with  a  frequency  tripled 
ND+J:  glass  laser  facility.  Optics  have  been  fabricated  which  will  allow 
line  focus  experiments  to  be  performed  with  this  facility  in  1982.  Design 
studies  on  soft  x-ray  cavity  reflectors  have  continued. 
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1.0  SUMMARY  OF  RESEARCH  OBJECTIVES 


During  the  January  1  through  December  31,  1981,  support  period 
we  have  concentrated  our  research  in  the  following  two  areas: 

1.  Extension  of  the  results  of  previous  studies  which  demonstrated 
the  existence  of  soft  x-ray  population  inversions  in  recombining, 
laser-produced  plasma  to  the  development  of  a  practical  amplifier 
device. 

2.  A  study  of  the  performance  which  would  realistically  be  expected 
to  be  obtained  from  multilayer  structures  suitable  for  use 

as  normal  incidence,  soft  x-ray  cavity  mirrors. 

Our  activities  during  this  period  have  led  to  several  significant  achieve¬ 
ments: 

1.  Demonstration  of  substantially  enhanced  population  inversion 
levels  in  point  focus,  step-target  geometry  due  to  the  use 
of  frequency  tripled  glass  laser  radiation  to  initially  heat 
the  plasma  medium. 

2.  Construction  of  cylindrical  focussing  plates  for  the  ultra¬ 
violet  (frequency  tripled)  laser  pulses  to  enable  line  focus 
experiments  to  begin  in  1982. 

3.  Successful  modelling  of  the  pc..  us  experimental  observations 

in  aluminum  plasmas  using  the  one-dimensional,  laser  plasma 

hydrodynamic  computer  code,  LILAC. 
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4.  Modelling  of  the  rol®  of  various  probable  fabrication  defects 
in  the  predicted  performance  of  multilayer  reflector  devices 
for  soft  x-ray  laser  cavities. 

An  elaboration  of  these  achievements  will  constitute  the  main  body 
of  this  report. 


2.0  STATUS  OF  THE  RESEARCH 


2.1  Soft  X-Ray  Amplifier  Development 

The  essential  idea  behind  our  approach  to  building  a  soft  x-ray 

laser  is  to  create  a  hot,  high  density,  highly  ionized  plasma  which  is 

converted  into  a  recombining  plasma  in  a  suitable  expansion  geometry. 

We  create  the  initial  plasma  conditions  by  focussing  a  high  intensity 

+3 

light  pulse  from  a  Nd  : glass  laser  onto  a  suitable  solid  target. 

During  the  past  18  months  a  significant  new  development  in  laser  tech¬ 
nology  has  occurred  at  LLE  which  has  greatly  enhanced  the  measured 
performance  of  such  recombining  plasmas  in  producing  inverted  popu¬ 
lations.  This  development  is  the  demonstration  and  now  routine  oper¬ 
ation  of  a  high  efficiency  frequency  tripling  system  for  high  peak  power 
+3 

Nd  rglass  lasers.  As  the  benefits  of  the  use  of  frequency  tripled, 
ultraviolet  radiation  in  these  experiments  became  clear,  we  modified  our 


experimental  facility  for  its  use.  In  this  section  we  will  review  the 


A 


frequency  tripling  concept  and  compare  the  use  of  ultraviolet  vs  infrared 
radiation  to  heat  high  density  plasmas.  We  will  then  directly  compare  / 


the  use  of  the  two  radiations  in  our  soft  x-ray  population  inversion 


experiment. 
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An  example  of  the  laser  target  irradiation  geometry  and  the  diagnostic 

set-up  we  use  is  given  in  Fig.  1.  In  a  typical  experiment  a  20-50 

+3 

‘Joule  0.5  nsec  FWHM  nd  :  glass  laser  pulse  is  focussed  onto  the  sur¬ 
face  of  a  moderate  Z  slab  target  such  as  aluminum  or  magnesium.  The 
target  geometry  is  arranged  to  allow  the  expanding,  hot  plasma  from  the 
surface  of  the  slab  to  be  partially  intercepted  by  a  foil  (shown  here 
as  magnesium). 

The  foil  is  postulated  to  provide  a  heat  sink  in  the  plasma  ex¬ 
pansion,*  leading  to  an  enhanced  rate  of  collisional  recombination  for 
the  highly  stripped  plasma  ions.  Such  recombination  processes  favor  the 

production  of  states  of  high  principal  quantum  number,  leading  to  pop- 
2 

ulation  inversion.  The  inversion  process  is  shown  schematically  in 
Fig.  2.  During  the  past  year  we  have  conducted  experiments  under  con¬ 
ditions  which  lead  to  a  substantially  increased  density  of  inverted 
population  on  4-*-3  transitions  in  helium-like  ions  of  the  slab  material 
(e.g.  aluminum  or  magnesium).  These  experiments,  described  in  the  next 
section,  greatly  increase  the  chances  for  a  successful  soft  x-ray  gain 
demonstration  in  the  next  year. 

The  frequency  tripling  scheme  has  been  described  in  a  series  of 

3  4  5  +1 

papers.  ’  *  The  linearly  polarized  output  of  a  Nd  :glass  laser  is 

incident  on  a  Type  II  KDP  crystal  which  is  angle-timed  for  optimum 

phase  matching.  The  input  polarization  is  set  at  35°  with  respect  to  the 

Ordinary  ray  direction  in  the  crystal  instead  of  the  usual  45®,  resulting 

.in  SH6  efficiencies  near  67%  for  a  broad  range  of  input  intensities. 
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The  light  emerging  from  the  doubling  crystal  thus  consists  of  a  1:1 
photon  mix  of  fundamental  and  second  harmonic  with  orthogonal  polarizations 
These  radiations  then  enter  a  Type  II  KDP  mixing  crystal  with  optimum 
mix  and  polarization  for  efficient  third  harmonic  generation.  The 
crystal  and  polarization  vector  orientations  are  shown  in  Fig.  3. 

Routine  overall  conversion  rates  to  third  harmonic  of  60%  or  more  are 
achieved  at  127  mm  aperture  in  the  system  used  on  our  glass  development 
laser  facility.  In  0.5  nsec  pulses,  between  40  and  50  Joules  of  0.351  ym 
radiation  are  delivered  to  targets,  resulting  in  on-target  intensities 
of  1015  W/cm2. 

A  systematic  study  of  many  aspects  of  laser  target  interaction 

physics  at  0.35  ym  has  been  conducted  over  the  past  year.  Many  results 

of  these  studies  are  published,6,7  or  in-press.  We  will  sunmarize  the 

relevant  features  of  ultraviolet  (third  harmonic)  irradiation  compared 

to  infrared  irradiation  of  targets  under  otherwise  similar  conditions. 

15  2 

For  on-target  intensities  in  the  vicinity  of  10  W/cm  we  find 

1.  The  UV  absorption  is  much  higher  than  the  IR  absorption 
(typically  75%  vs  30%). 

2.  The  UV  absorption  is  classical  (collisional )  as  opposed  to 
resonance  absorption  in  the  IR.  There  is  no  suprathermal 
electron  distribution  or  suprathermal  energy  transport  with  UV 
irradiation. 
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3.  The  UV  absorption  occurs  at  higher  electron  density  1022  cm-'*) 
than  the  IR  absorption  1021  cm-3). 

4.  The  plasma  coronal  electron  temperature  is  the  same  or  slightly 
higher  with  UV  irradiation  (1-2  keV). 

5.  Both  soft  x-ray  line  and  continuum  emission  increases  by  an 
order  of  magnitude  or  more  with  UV  irradiation  compared  Joule 
for  Joule  with  the  IR. 

These  observations  suggested  that  frequency  tripled  Nd  :glass 
laser  radiation  would  be  advantageous  for  use  in  our  population  inversion 
experiments.  In  particular,  it  appeared  that  we  could  «ary  one  of  our 
initial  conditions  -  the  plasma  particle  density  -  while  leaving  the 
initial  temperature  and  the  expansion  geometry  unchanged.  The  absence 
of  hot  electrom  production  would  also  help  to  carify  the  role  of  the 
plate  we  employ  in  the  expanding  plasma.  If  the  role  of  the  plate  is 
strictly  geometricaly  or  thermal  in  nature  then  it  is  possible  that  our 
density  of  inverted  population  might  be  increased. 

To  implement  the  conversion  to  frequency  tripled  light  in  our 
experiment,  our  laser  beam  transport  optics  were  recoated  for  0.351  urn, 

the  vacuum  window  to  our  experimental  chamber  was  replaced  by  a  coated, 
fused  silica  window,  and  our  focussing  lens  was  replaced  by  a  single 
element  f/ 12  fused  silica  lens.  Alignment  of  the  UV  system  is  accomplished 

with  the  aid  of  a  CW  argon  ion  laser  operating  at  0.355  ym.  This  CW 
'beam  is  Inserted  into  the  beam  transport  system  at  the  point  where  the 
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IR  laser  enters  the  frequency  conversion  system.  A  schematic  diagram  of 
the  frequency  tripling  system,  the  alignment  system  and  the  beam  diagnostics 
is  given  in  Fig.  4. 

We  now  present  some  comparative  results  on  step-target  experiments 
using  IR  and  UV  irradiation  under  otherwise  similar  conditions.  In 
these  experiments  a  target  geometry  similar  to  that  shown  in  Fig.  1  was 
used  except  that  the  foil  was  arranged  in  the  form  of  two  half-planes, 
i.e.  a  slit.  The  spatially  resolving  crystal  spectrograph,  set  to 
record  transitions  to  the  ground  state  of  helium-like  and  hydrogen-like 
target  ions,  was  positioned  to  view  along  the  open  axis  of  the  slit. 

In  Fig.  5  we  present  two  spatially  resolved  spectra  from  aluminum 
targets  heated  by  50  Joules,  600  psec  FWHM  pulses  of  1.054  ym  laser 
radiation,  recorded  using  a  flat  TAP  crystal  spectrograph  with  a  spatially 
resolving  slit.  The  equivalent  source  space  dimensions  for  our  geometry 
are  shown  on  the  figure.  A  lead  foil  slit  was  used  in  the  target 
producing  the  left  hand  spectrum.  This  foil  was  convenient  in  that  it 
produces  no  strong  emission  lines  in  this  part  of  the  spectrum.  A 
similar  expansion  spectrum  is  observed  in  the  right  hand  spectrum, 
obtained  from  an  aluminum  target  fitted  with  a  thick  magnesium  slit  in 
the  form  of  a  channel.  Because  of  the  close  proximity  of  many  of  the 
magnesium  and  aluminum  lines  is  it  difficult  to  perform  line  ratio 
analysis  for  temperature  and  density  profile  reconstruction.  However, 

'the  use  of  lead  foil  targets  consistently  produced  severe  damage  on  the 
laser  focussing  optics  due  to  production  of  high  velocity  lead  debris 
along  the  plasma  expansion  direction. 
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In  Fig.  6  we  present  two  spatially  resolved  spectra  from  aluminum 
targets  heated  by  20  Joule,  490  psec  FWHM  pulses  of  0.351  urn  laser 
radiation,  recorded  on  the  same  crystal  spectrograph  as  used  to  record 
Fig.  5  except  that  an  additional  0.75  y  aluminum  foil  filter  reduces 
the  greatly  increased  x-ray  emission  intensities  with  the  UV  irradiation. 
The  equivalent  source  space  dimensions  are  shown  on  the  figure.  The 
left  hand  spectrum  was  produced  by  a  target  with  a  titanium  foil  slit 
while  the  right  hand  spectrum  was  produced  by  a  target  with  a  magnesium 
foil  slit.  A  slight  misalignment  of  the  spatially  resolving  crystal 
spectrograph  axis  with  respect  to  the  open  axis  of  the  target  slit  is 
evident  in  this  figure.  The  two  spectra  show  very  similar  expansion 
characteristics,  indicating  that  the  performance  of  these  targets  in  the 
UV  is  essentially  independent  of  the  heat  sink  material  just  as  has  been 
observed  for  IR  irradiation. 

The  spectra  shown  in  Figs.  5  and  6  are  presented  on  nearly  identical 
spatial  scales  and  recorded  optical  densities.  It  is  immediately  striking 
that  the  length  of  expanding  plasma  column  giving  rise  to  intense  radiation 
from  the  n  ■*  1  transitions  is  much  greater  with  UV  target  irradiation. 

This  is  consistent  with  the  coupling  of  UV  radiation  and  hot  plasma 
formation  at  much  higher  particle  density  than  with  IR.  Densitometric 
scans  of  the  spectra  in  Fig.  6  show  a  clearly  defined  inversion  in  the 
intensity  of  the  4  1  over  the  3  -*>  1  transitions  in  the  helium-like 
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SHOT  NUMBERS  2789  AND  2788 


aluminum  ion  beginning  between  600  and  900  um  from  the  initial  target 
surface,  roughly  2-3  times  farther  out  than  is  seen  in  the  IR  experiments. 

Based  on  the  absolute  intensity  calibration  of  our  spectrograph  we  calculate 
the  density  of  inverted  population  on  the  4  -+  3  transition  would  be  3-5 
times  higher  in  these  UV  experiments  than  in  any  of  our  IR  measurements. 

As  a  result  of  our  decision  to  convert  our  experimental  system  to 
ultraviolet  irradiation,  we  postponed  our  procurement  of  a  cylindrical 
focussing  lens  for  line  focus  experiments.  The  improved  performance  of 
the  recombining  plasmas  with  ultraviolet  irradiation  has  motivated  us 
to  redesign  our  cylindrical  focussing  system  for  the  ultraviolet.  The 
f/12  spherical  focussing  lens  has  a  focal  length  of  approximately  1660  mm. 

This  is  sufficiently  long  to  permit  the  use  of  additional  weak  cylindrical 
plates  to  obtain  a  line  focus.  By  employing  two  such  plates  we  may  obtain 
continuous  control  over  the  aspect  ratio  of  the  line  focus  by  rotation  of 
their  respective  cylinder  axes.  A  pair  of  fused  silica  plates  has  been 
fabricated  for  this  purpose.  These  plates  are  currently  deployed  in 
preliminary  line  focus  experiments. 

Because  of  the  essentially  classical  behavior  of  intense,  focussed, 
ultraviolet  laser  pulses  in  plasma  media,  we  felt  that  for  the  first  time 
we  might  hope  to  model  the  experiments  on  production  of  population  inversion 
described  above.  Accordingly,  in  our  first  such  efforts  we  sought  to  adapt 
a  highly  developed  numerical  plasma  hydrodynamic  code,  LILAC,  to  the 
description  of  the  observed,  time  integrated,  spatially  resolved  emission 
spefctrum  from  the  plasma  ions.  LILAC  calculates  inverse  bremstrahlung 
absorption  for  an  electromagnetic  wave  propagation  through  a  one  dimensional 
(i.e.  spherically  symnetric)  plasma.  The  hydrodynamic  equations  of  continuity. 
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and  energy  and  momentum  conservation  are  then  solved  at  each  point  during 
each  small,  adjustable  time  interval.  The  temperatures  of  electron  and 
ion  populations  are  calculated  separately,  and  energy  exchange  between 
populations  is  provided  for  in  each  computational  cell.  The  various  coronal 
species  are  assumed  to  behave  in  an  ideal  gas  fashion.  The  degree  of 
ionization  is  calculated  at  each  step  using  the  model  appropriate  to 
the  plasma  conditions. 

To  enable  a  one-dimensional  code  to  accurately  replicate  observations 
made  on  flat  target  experiments,  it  has  been  found  effective  to  assume 
an  irradiation/target  geometry  shown  in  Fig.  7.  If  the  diameter  of  the 
laser  focal  spot  on  the  flat  target  is  $  then  a  spherical  target  radium 
of  2« t  is  chosen  in  the  calculation.  The  expansion  of  the  plasma  proceeds 
radially  from  the  surface  in  a  pie-shaped  region. 

The  emission  spectrum  is  obtained  using  a  rate  equation  package 
which  is  calculated  as  a  post- processing  routine  on  LILAC.  The  excitation 
and  recombinaton  coefficients  for  the  helium-like  and  hydrogen-like 
species  are  adapted  from  McWhirter  and  Hearn.  Only  ionization  ground 
state  populations  are  calculated  for  lower  stages  of  ionization.  The 
calculations  are  suppressed  at  all  points  for  which  the  electron  temper¬ 
ature  drops  below  60ev  due  to  the  lack  of  good  data  available  for  the 
species  and  levels  of  interest  at  such  energies. 

To  replicate  the  observations  of  our  spatially  resolving  spectrograph, 
the  expanding  plasma  is  divided  into  a  series  of  laminar  layers,  10  ym 
thick,  oriented  in  the  direction  of  observation.  The  radiation  from 
each  laminar  segment  (assumed  to  be  distinguishable  in  our  instrument) 

Is  calculated,  taking  into  account  the  opacity  of  each  ground  state 
transition  using  the  model  of  Elton.  Each  simulation  yields  a  time-  and 
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space- resolved  spectral  emission  distribution  as  well  as  a  time-integrated, 
space- resolved  distribution. 

To  simulate  the  effect  of  the  foil  structure  in  the  step- target 
experiments,  we  add  a  spatially  dependent  energy  loss  term  to  the  hydro- 
dynamic  equations  in  the  code,  in  thich  the  cooling  tine,  tc,  is  an  adjustable 
parameter.  In  this  way  the  results  could  be  compared  to  a  particular 
physical  model  of  the  transfer  of  energy  from  the  plasma  to  the  cooling 
plates  or  foils. 

In  Fig.  8  we  present  a  calculation  of  the  time- integrated,  spatially 
resolved  emission  from  the  4  -*•  1  and  3  1  transitions  in  helium-like 

aluminum.  In  this  simulation,  a  peak  irradiation  intensity  of  5  x  10lu 
W/cm2  at  0.351  pm  was  assumed  with  a  focal  spot  diameter  of  100  pm.  A  cooling 
plate  100  pm  thick  was  located  200  pm  from  the  initial  target  surface.  The 
laser  pulsewidth  was  500  psec  and  assumed  to  be  gaussian. 

A  "normal"  line  intensity  distribution  is  exhibited  close  to  the  initial 
target  surface.  At  a  distance  of  approximately  700  pm  from  the  target  the 
two  intensities  have  become  equal  and  beyond  that  point  they  are  inverted. 

This  behavior  is  in  very  good  agreement  with  the  experimental  observations 
presented  earlier.  To  obtain  this  result  a  cooling  time,  tc,  of  30  psec 
was  used  in  the  calculation.  It  appears  that  the  exact  value  chosen  for 
t  is  not  critical.  For  t  =  10  psec,  the  increase  in  line  emission  in 

C  t 

the  vicinity  of  the  cooling  plate  appears  to  be  somewhat  larger  than 
observed.  For  t  >  100  psec  the  lower  temperature  drop  substantially  reduces 

w 

the  emission  from  this  region. 

The  initial  success  of  this  modelling  effort  is  very  encouraging. 

During  the  current  support  period  we  plan  to  continue  this  effort  to  see 
If  further  refinement  and  optimization  of  the  experimental  conditions  can. 
be  verified. 
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STEP  TARGET  SIMULATION  iL I  LAC 


DISTANCE  PROM  THE  TARGET  (MM) 


2.2  Soft  X-ray  Cavity  Development 


Since  we  began  our  study  of  periodic  structures  which  could  be 
-suitable  for  use  as  reflecting  elements  in  a  soft  x-ray  laser  cavity, 
considerable  activity  has  been  initiated  in  several  laboratories  in 

O 

the  development  of  such  elements.  Substantial  variation  in  the 
technical  approach  to  fabrication  of  such  elements  currently  exists 

Q 

and  at  least  one  attempt  at  commercial  manufacture  may  soon  be  realized. 
During  the  past  year  we  have  generalized  our  analysis  of  the  performance 
which  could  be  expected  from  such  elements  to  include  a  treatment  of 
the  effects  of  various  structure  defects  and  fabrication  errors.  Interim 
results  of  this  work,  which  is  continuing  at  present,  are  presented  in 
the  attached  publication  reprints  and  present  a  fair  summary  of  the  work 
performed  during  the  support  period. 


i 
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ABSTRACT 

We  treat  o  variety  of  problems  in  the  theory  of  soft  x-roy  multilayer  mirrors,  A 
characteristic  matrix  solution  to  Maxwell's  equations  is  presented  that  applies  to  both 
periodic  and  non-periodic  reflectors  whose  layers  con  possess  arbitrary  index  gradients. 

Procedures  ire  derived  to  maximize  multilayer  reflectivity  in  the  kilovolt  and  sub-kilovolt 

regime.  A  refractive  correction  to  Brogg's  law  is  derived  thot  includes  the  effect  of  < 

absorption  as  well  as  the  effect  of  dispersion.  Multilayer  reflectivity  in  the  presence  of 

random  thickness  errors  is  treated  analytically.  An  analytic  treatment  of  different  kinds 

of  interfacial  roughness  is  described.  The  reflecting  properties  of  the  multilayers  moy 

contain  qualitative  signatures  that  ere  characteristic  of  the  different  kinds  of  roughness. 

The  effect  of  interlayer  diffusion  is  discussed. 


I.  Characteristic  Matrix  Analysis  of  X-Ray 
Reflectors 


Our  analysis  of  x-roy  multilayers  proceeds  from  a 
solution  to  Maxwell's  material  equations,  given  that 
the  multilaver  structure  has  a  spatially  varying 
complex  dielectric  constant  ♦.  In  this  respect  our 
analysis  follows  the  Ewald-von  Laue  dynamical  theory  of 
crystal  diffraction.  In  the  dynamical  theory  it  is 
assumed  that  <  is  three-dimensionally  periodic,  as  in  a 
perfect  crystal.  We  assume  that  •  is  a  function  only 
of  the  coordinate  z  normal  to  the  multilaver  surface;  < 
need  not  be  rigorously  periodic. 

In  this  case  it  is  known  that  the  wave  equation 
can  be  separated  into  ordinary  differential  equations, 
whose  solutions  can  be  put  into  characteristic  matrix 
form.1  The  characteristic  matrix  solutions  for 
homogeneous  layers  are  commonly  used  in  optical 
multilayer  calculations. 

In  the  visible,  these  differential  equations  must 
be  solved  numerically  in  the  case  of  a  eneral 
structural  profile,  but  for  the  x-ray  regime  we  have 
solved  them  analytically  under  the  assumption  that 
|  i  -I  |  «  I  .  Essentially,  we  treat  the  material  medium 
as  giving  rise  to  a  perturbation  in  the  vacuum  fields. 

We  now  present  this  solution. 

Let  the  plane  of  incidence  be  the  y-z  plane. 


Define: 

S  Polarization  P  Polarization 


E  ;E* 

H  =  H  /cos8 
y 


E  5  -§  /cos9 
y 

H  s  A 


(1) 


The  characteristic  matrix  solution  for  the  Kth 
cell  of  the  multilayer  shown  at  right  is: 


•cos*c-pt 
4  (t«+-yx) 


-i  (ta-yx) 

«*  #**p* 


(2) 


where: 

A**' cos f 

<«*>* 


y* ■  JZ.sk  • /d* A<z>eoi[_iL  costs] 
*  ctU  V  *  I 

■  -T  *  *  *•»  •*  (t-  «"  •*) 
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Tig.  1.  Decomposition  of  reflector  into  wit  cello. 


We  note  that  the  structure  as  determined  by 
Alt)  ■  |/2(«UH)  need  not  correspond  to  o  succession  of 
homogeneous  layers,  but  can  have  an  arbitrary  variation 
along  z. 

If  ph  is  the  amplitude  reflectivity  of  the  stack  of 
cells  K-l,  K-2,  ...  I,  them 


l  +  pa  Htm  t-pp  W 


From  these  we  can  derive  a  difference  equation  . 
that  propagates  the  amplitude  reflectance  from  wit 
cell  to  unit  ceil.  This  is  analogous  to  the  amplitude 
recursion  formula  that  is  used  to  propagate  the  « 
reflected  amplitude  from  single  layer  to  single  layer 
in  optical  multilayer  calculation*.  The  equation  fan 


w  soc  #  /  di  Ad) 

*  ort 


tx  ■  ue  e«  • 


(3) 


This  equation  is  particularly  useful  in  analysing 
non-periodic  multilayers. 


(Hoe  wt  fen lee  (tw  optics  ceevmtiee  eherc  f  is  mestund  to  the  ootmoL) 
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IL  Equivalent  Parameter  Analysis  of  Periodic 
Reflectors 


where 

This  is  olso  the  reflectivity  of  at  x-roy 
rrultiloyer  having  J  cells. 


Fig.  2.  Structure  of  tetit  cell  in  e  bileyer  reflector. 


A  multilayer  with«a  bilayer  structure  consists  of 
alternating  homogeneous  layers  of  two  different 
materials.  If  the  thicknesses  of  the  different  layers 
of  each  material  ere  equal,  the  multilayer  is  periodic. 

By  placing  the  cell  boundaries  at  the  midpoints  of  the 
L  layers  the  cell  structure  con  be  made  centrosymmetric 
as  shown  in  the  diagram  above. 

Then  the  parameter  p  is  zero,  and  the 
characteristic  matrix  solution  becomes: 


(- 


(6) 


where: 


q 


'-(t) 

■(f)' 


2dV  (Ah  •  Ai)  sin** 


[«A I  +  A-  (Ah  -  All] 


(7) 


This  has  the  form  of  the  well-known  characteristic 
matrix  solution  for  a  single  homogeneous  layer  having 
«n  index  of  refraction  a.  and  phase  thickness  A : 1 


‘V 

J 

\ 

'HO.  mb  a 

Hi /n.)  tin  A 

•«A 


(8) 


III.  Optimization  of  Reflectivity 

A  computer  progrom  con  be  used  to  determine  the 
optimum  layer  thicknesses  in  an  x-roy  reflector  vio 
standard  optimization  algorithms.  A  study  of  such 
optimized  designs  shows  thot  if  the  radiation  energv  is 
above  about  100  ev,  a  bilayer  reflector  that  is 
periodic  can  achieve  neorly  the  maximum  reflectivity 
thot  is  passible  from  given  layer  materials. 

The  results  of  the  computer  study  ogree  with  a 
theoretical  aolysis  we  hove  mode.  It  can  be  shown 
analytically  thot  if  J  is  large,  the  reflectivity  of  an 
optimized  periodic  reflector  is  an  extremum  with 
respect  to  an  arbitrary  perturbation  in  structure. 

(The  structures  of  the  different  cells  need  not  be 
perturbed  equally).  This  onolysis  neglects  smoll  terms 
proportional  to  the  squore  of  the  difference  of  the 
index  of  refroction  from  one. 


NUMBER  OF  LAYER  PAIRS  J 

1  q (j) ei me putettc d— H«» Hpttmiw w—ctiz  pjinniitiiHW" 
1.  q  (J)  ¥  tni$m  o»«nm»d  «r  wee  J  (2qr«oHiH»»i 

1.  q  (J)  •!  •  perwdtc  aniqn  oqWniMe  M  J  —  •>  (I  e»qw«t  It  treeeowi 


Fig.  J.  Comparison  of  optimization  schemes. 


To  motch  this  matrix  to  the  characteristic  matrix  for 
the  unit  cell  of  an  x-roy  reflector,  eq.  (6). we  must  moke  the 
assignments  A  "  «r  r - yl  ax)  <>->)  ,  where 

t  m  .  The  two  motricies  will  be  equal  to  within 
first  order  in  A  wtd  A. 

This  matching  of  motricies  is  ai  exomple  of  whot 
in  thin,  film  optics  is  known  as  the  method  of 
equivoifctt  parameters.1  There  is  o  straightforward 
analogy  between  the  proparties  of  the  single  equivalent 
layer  and  the  single  cell  of  x-roy  reflector. 

A  reflector  with  J  cells  corresponds  to  a  stack  of 
J  equivalent  layers,  which  is,  in  effect,  o  single 
layer  of  phase  thickness  if.. 

By  letting  £,■!  +  *,  and  N,  ■  I  -  e,  ,  we  can 
obtain  the  well-known  egression  far  the  reflectivity 
of  a  single  layer: 


„  _  Ijjk-UAMeejlg 
0m  2  •  ks.  ♦  I/A)  tee  UR 


(9) 


For  this  reoson  we  now  concentrate  on  the  periodic 
case.  When  J  is  lorge,  the  reflectivity  of  a  periodic 
multiloyer  with  a  centrosvmmetric  structure  is  give-  bv 

p  ■  -  —  where  8  at/r'  -  V  U0' 

y 

Let  o  dot  represent  a  derivative  with  respect  to 
some  Structural  parameter.  Using  q« -(•+!)  >  one  con 
obtoin  o  general  optimization  condition  on  the 
intensity  reflectivity 

».»  *,(£).»  m(^)-o  up 

A  similar  but  slightly  more  complicated  result  con 
be  obtoined  if  prO. 

We  show  below  in  sec.  IV,  eq.  16-17,  thot  ot  the 
wavelength  of  peak  reflectivity,  I  is  pure  imoginorv. 
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Rc  0/«>  -  0 


The  condition  for  maximizing  the  peak  reflectivity 
is  therefore: 

or 

ln(i-t  —  )»0  (12) 

y 


Rc  <«>  -  0 


(17) 


As  an  application  of  this  result,  we  consider  the 
problem  of  optirriring  the  H  to  L  thickness  ratio  in  a 
multilayer  with  a  bilaver  structure. 

In  order  to  do  so  we  apply  our  optimization 
condition  to  the  parameter 

'ST£ 

and  obtain  a  condition  first  obtained  by  Vinogradov  and 
leldovich  for  the  case  of  reflection  at  normal 
incidence:’ 


tan  &.»».+  jrltnIAJ  (t4) 

la  (Ax  -  Al) 


More  generally,  we  consider  mth  order  diffraction, 
jid  let  fit,  “  a*  (da  dN+di.).  We  obtain  the  optimization 
condition: 


uik'fct 


imrlmCAi.) 

I  m  (An  -  At.) 


(15) 


In  the  plot  below,  and  in  the  other  numerical 
examples  of  this  report,  ’ve  have  used  preliminary 
values  for  atomic  scattering  factors  from  a  forthcoming 
compilation.*  These  have  generously  been  provided  to 
us  by  B.Menke. 


The  some  condition  is  obtained  in  the  general  case 
when  p*0. 

The  requirement  thot  Rc  (*)  “  0  is  the  same  as 
requiring  thot  the  real  part  of  the  equivalent  phase 
thickness  A  equal  w,  which  is  different  from  the  usual 
requirement  in  the  absence  of  absorption  thot  the  reol 
part  of  the  optical  phase  thickness  equal  ».  The 
difference  is  a  consequence  of  the  presence  of  multiply 
reflected  beams.  In  an  absorbing  structure,  these 
beams  undergo  phase  changes  upon  reflection,  causing 
phase  change  in  the  overall  oscillation  across  eoch 
cell.  r*— 

Since  4  ■  Vr  -  y!  ,  we  con  write  our  requirement 
that  Re  (4)  *  0  as  Im  (41)  -  0  ,  which  can  be  shown  using 
the  Schwarz  inequality  to  automatically  imply 
Re<4!)  <  0  .  The  condition  In  (i‘)  •  0  can  be  rrxnipuloted 

to  give: 

(-•  <>»> 


where  primes  aid  double-primes  denote  real  aid 
imaginary  pcrts.  This  result  was  first  discovered  by 
F.  Miller’  in  the  context  of  crystal  diffraction,  rod 
independently  by  us  in  the  context  of  multilayer 
reflection.  P.  Lee,  of  LASL  is  a  collaborator  with 
us  in  this  work. 

In  the  usual  formula  for  the  dispersion 
correction,  the  term  in  square  brackets  is  missing.* 
The  first  term  in  parentheses,  it' ,  can  be  considered 
to  represent  dispersion,  the  term  in  square  brackets 
absorption. 

The  relative  magnitude  of  the  two  effects  are 
compared  in  the  plot  below.  The  ratio  plotted  is 


Ratio  ■ 


8i-8, 

e„-e2 


(19) 


where  ft*  arcos  (A/2d),  ft  is  the  true  Bragg  angle  os 
corrected  for  absorption  and  dispersion,  and  ft  is  the 
Bragg  angle  as  corrected  for  dispersion  only. 

It  can  be  shown  thot  the  absorption  correction  is 
always  smaller  than  the  dispersion  correction,  but  as 
the  plot  shows  there  are  cases  in  which  the  two  are 
comparable.  The  absorption  correction  must  therefore 
be  included  whenever  on  accurate  matching  of 
wavelength,  angle,  and  2d-spacing  is  required. 


Fig.  4.  1  'sftectivity  as  a  function  of  thickness  ratio. 


IV.  The  Absorption  Correction  to  Bragg's  Law 


It  is  well-known  thot  the  Bragg  condition  for 
crystal  diffraction.  Id  cos  •  *  X  ,  must  be  corrected  for 
the  dispersive  effects  of  the  crystalline  medium 
(following  the  optics  convention,  ft  is  the  angle  of 
incidence  to  the  normal).  It  is  less  well-known 
however,  thot  Bragg's  low  must  be  corrected  for 
absorption  as  well  as  far  dispersion. 

In  our  forrrfalism,  the  parameter  ft  represents  the 
detuning  of  a  multilayer  from  the  simple  Bragg  angle 
ft  a  areas  (4  Id),  through  the  relation: 

♦  *  (eos9»  -  eos8  )  (16) 

To  find  the  exact  angle  of  maximum  reflectivity  we 
optimize  the  structure  with  respect  to  the  parameter  ft. 
Then  using  the  optimization  condition  introduced  above 
<eq.(ll)),  and  setting  y  -  0,  i  -  I,  we  get: 


Fig.  5.  Comparison  of  abrarptlon  correction  to  dispersion 
correction. 
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V.  The  Effect  of  Random  Thickness  Errors  On 
Multilayer  Reflectivity 


<  A  >“ 


I  ♦  »Ei 
I  •  »£> 


•  20 
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Rmdom  thickness  errors  in  tn  x-ray  multilayer  con 
awes  o  severe  degradation  in  reflectivity.  Even  when 
the  errors  ere  small  compared  to  the  multilayer's 
Idk^acing,  they  con  cause  a  significant  cumulative 
defeating  if  not  compensated  for  during  the  fabrication. 

Previous  theoretical  molysis  of  the  effect  of 
tsrtdom  thickness  errors  is  largely  limited  to  the  work 
of  Shellan  an  optical  Brogg  reflectors.7  She II an 
treats  thickness  errors  as  a  small  perturbation  in  the 
Structure  of  dielectric  reflectors  having  a  small 
coupling  constant  per  ceil.  X-roy  reflectors  have  o 
■mall  coupling  constant,  but  it  may  not  be  adequate  to 
Consider  the  errors  in  x-ray  reflectors  to  be  small 
pvturtoatians  that  cause  only  a  small  degradation  in 
reflectivity.  In  the  soft  x-roy  regime  It  is  also 
important  to  include  absorption,  which  can  be  the  major 
factor  determining  the  reflectivity. 

We  have  met  these  requirements  using  an  analysis 
boMd  an  the  amplitude  recursion  equation  introduced 
above  (aq.  (5))l  This  equation  is: 

p  v  1  p-(iy-p)-( *y*p)  p*  (20) 

K*|  K  K 

We  treat  Ok  as  a  random  variable  with  mean 

<d»>  ■e-(Je<4,>/X)cos»  (21) 

and  known  variance  <  A*’  >.  In  o  bilayer  reflector, 

<  AO1  > ““p  «■*’ * ( <  AL*  >  +  <  AH1  >  )  (22) 

If  the  errors  in  the  L  and  H  layers  are  independent. 

We  make  the  problem  analytically  tractable  by 
first  dividing  ox  into  a  deterministic  aid  □ 
nan-deterministic  pat,  Ox  ■  <  a*  >  +  A  ,  and  then 
making  use  of  the  inequality: 

<l*l*>  <  <|e*lI>  <  I  (23) 

The  r tgpnt  inequality  is  fairly  strong  in  the  soft 
x-ray  regime,  where  absorption  is  fairly  large.  The 
rifftt  inequality  is  olso  p«ficularly  strong  when 

<  Ad’  >  Kiel,  since  the  reflectivity  is  then 
severely  degraded  by  the  thickness  errors. 

On  the  other  hand,  the  left  bond  inequality 
becomes  very  strong  when  <  Ad'  >  <  (  o  ( ,  since 

<  Ox  >  ■>  Ox  when  the  errors  are  smal  I. 

For  these  reasons  l)xl  is  generally  ■*  |;  in  our 

aialysis  we  therefore  neglect  terms  that  are  cubic  or 
hitter  in  pk-  Our  approximation  is  accurate  in  the 
limits  <  Ad1  >  ►  I  el  and  <Ad*>dlo!,  but  is  olso 
fairly  accurate  in  the  intermediate  region. 

In  a  simila  way  we  treat  the  K -dependence  of  Ox 
md  ps  in  at  approximate  way  that  is  correct  in  the 
limits  <Ad’>  d  I  el  i<Ad*>  H el.  X  *  I  el'1,  or  when 
the  absorption  is  sufficiently  lags  that  |  «|  <  I. 

The  finol  result  fa  o  reflector  with  J-]  cells 
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The  accompanying  figure  shows  < X,  >  plotted  os  o 
function  of  RMS  thicness  error  in  o  paticular  example. 

We  have  tested  the  oecurocy  of  the  analysis  by 
moking  o  Monte  Carlo  simulation.  A  computer  progrom  is 
used  to  generate  a  large  number  of  randomly  perturbed 
multilayer  stocks,  and  to  compute  their  mean 
reflectivity.  Sample  results  are  shown  in  the  plot. 

The  plot  olso  shows  o  computation  of  AA  k  as  o  function 
of  RMS  thickness  error. 

There  ore  no  free  parameters  in  the  fit  of  the 
analytic  egression  to  the  Monte  Carlo  results. 


Fig.  6.  Effect  of  random  tnicknets  error*  on  multileyer 
reflertivity. 


VI.  The  Effect  of  Interfocial  Roughness  On 
Multilover  Reflectivity 


Eastman  has  developed  numerical  methods  to  treat 
the  effect  of  loyer  roughness  in  opticol  multilayers. • 

We  hove  used  similar  physicol  assumptions  to  treot 
eerloin  kinds  of  roughness  in  x-roy  multilovers 
malytically.  Our  model  assumes,  in  effect,  that  the 
near-field  reflected  beom  above  tnv  point  on  the 
multiloyer's  surfoce  con  in  principle  be  calculated  bv 
inserting  locol  volues  for  the  loyer  properties  into  an 
algorithm  that  computes  the  reflectivity  of  multilayers 
with  plonor  interfaces.  Such  a  scalar  model  of  the 
roughness  requires  that  the  tronsverse  autocorrelation 
length  of  the  roughness  be  forge  compared  to  the  loyer 
thicknesses.  A  more  deloiled  discussion  of  the 
physicol  assumptions  of  the  scolor  model  is  given  in 
ref.  B. 

Ref.  ?  discusses  the  limitations  of  scolor 
scattering  fheorv  in  comparison  with  more  rigorous 
tfwories.  In  general  the  scolor  theory  is  oest  at 
predicting  fotol  specular  and  diffuse  reflectivities 
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,jxd  at  predicting  the  angular  distribution  of  the 
diffusely  reflected  beam  at  angles  close  to  the 
specular  beam;  it  cannot  predict  polarization  effects. 

One  kind  of  multilayer  rouqhness  that  has  been 
treated  with  the  scolor  theory  is  that  which  Eostmon 
calls  "identical  films'",  in  which  all  layers  are 
considered  to  reproduce  a  common  roughness  profile 
Uenerally  that  of  the  substrate).  The  autocorrelation 
enqth  of  the  roughness  in  the  longitudinol  direction 
■  s  therefore  very  large.  Spiller  ef  at.  and  Haelbich 
et  al.  hove  discussed  the  effect  of  roughness  on  x-ray 
mult. lover  reflection  in  terms  of  the  same  expression 
as  results  from  the  identical  film  model.  1  **l‘ 

We  have  modeled  two  kinds  of  non-identicol 
-?i<jhness  In  which  the  longitudinal  autocorrelation 
Vnqth  is  very  small,  so  that  the  roughness  profiles  of 
me  different  lovers  are  uncorrelafed.  We  refer  to  the 
two  as  "rauqhening  films'"  and  "'smoothening  films'". 

In  the  case  of  roughening  films,  we  assume  that 
the  errors  in  the  local  layer  thicknesses  above  each 
point  on  the  surface  cause  a  cumulative  dephasing,  so 
•hat  the  absolute  rouqhness  of  the  top  layer  increases 
m  a  random  walk  fashion  as  more  layers  are  added.  One 
mcv  consider  the  formation  of  these  films  to  be  such 
that  the  granularity  In  each  layer  is  independently 
acded  to  a  baseline  of  roughness  established  by 
preceed.nq  layers. 

Jnder  the  assumptions  of  the  scalar  model,  the 
near-field  amplitude  is  given  by  the  equation  for 
one-dimensional  thickness  errors  discussed  above.  The 
far  f  eld  coherent  reflectance  is  obtoined  by 

evaluating: 


coh 


5  I  <<5f  ar  field  *  ^ 


J-l 

2 x) 

=  ]  Earn  er  fran$form(  <0  jq(  x  )  e  ‘  K 


-tere  x  is  a  coordinate  along  the  surface,  and  pio<*)  is 
•ne  necr-f  eld  amplitude  as  measured  at  the  upper 
Surface  of  the  multilayer.  This  upper  surface  is 
rough,  so  the  factor  cxp(2iXAe«)  must  be  used  to 
propagate  e.ofs)  to  a  mean  plane,  where  the  far-field 
amp.ifjde  can  properly  be  evaluated  via  the  Fourier 
transform. 

The  argument  of  the  transform  Is  independent  of  x, 
so  R.w  will  be  a  delta-function  of  the  angle  of 
reflection.  can  therefore  be  identified  as  the 

specular  beam.  Our  expression  for  R—  in  the  soft  x-ray 
regime  s  given  below  (eq.',27)). 

A  diffuse  beam  is  also  present.  The  total 
intensity  of  the  diffuse  and  specular  beams  is 
determined  by  the  total  absorption,  which  is  the  same 
as  in  the  case  of  JO  thickness  errors  discussed  above. 

We  have  also  developed  an  analytic  model  for  fhe 
kind  of  roughness  we  call  "smoothening  films”.  Such 
films  may  be  considered  to  have  a  leveling  nature 
during  some  stage  of  formation,  but  to  nonetheless 
possess  <F>  intrinsic  roughness  after  formation  is 
complete.  We  therefore  assume  that  an  error  in  the 
local  thickness  that  a  lover  has  ot  some  position  on 
the  reflector  will  be  compensated  for  in  the  thickness 
of  the  next  layer  deposited.  The  roughness  height 
necessary  to  cause  o  given  drop  in  reflectivity  is 
therefore  much  larger  than  in  the  cose  of  roughening 
films,  because  the  thickness  errors  do  not  accumulate 
(see  plot). 

It  turns  out  thot  the  smoothening  property  of  this 
type  of  roughness  causes  fhe  intensities  of  the  diffuse 
and  specular  beams  fo  become  equal  only  of  a  level  of 
roughness  where  the  total  reflectivity  has  been 
decreased  quite  substantially  (via  an  increase  in 
absorption).  In  contrast,  with  roughening  films,  fhe 
two  become  equal  at  o  roughness  level  where  the  total 
reflectivity  is  only  slightly  decreosed.  With 
identical  films,  the  total  reflectivity  is  unaffected 
by  the  magnitude  of  fhe  roughness. 

With  films  of  both  the  roughening  and  smoothening 
types,  the  total  absorption  reaches  a  steady-state 


level  as  more  avt  more  layers  are  added.  However  In 
the  cose  of  roughening  films  the  proportion  of  the 
reflected  light  in  the  specular  beam  steodity 
decreases,  since  the  upper  surfaces  get  steadily 
rougher. 

The  acceptance  ingle  of  multilayers  with 
smoothening  films  is  not  greotly  influenced  by  the  RMS 
magnitude  of  the  roughness.  In  the  case  of  roughening 
films,  the  acceptance  angle  is  increased  in  somewhat 
the  same  way  as  would  be  caused  by  at  increase  in  the 
layer  bulk  absorption  constants. 


Fig,  7.  Effect  of  roughness  -  reflection  into  specular  beam. 


Our  egressions  for  specular  reflection  in  the 
presence  of  roughness  ossume  thot  the  absorption  hos 
reached  its  steady-state  value.  Our  expression  for  the 
case  of  smoothening  films  assumes  a  bilayer  structure 
and  Coussicr  statistics  for  the  roughness. 

Our  egression  for  roughening  films  is: 

2  (F,  - G  ,-,) 

i#  ♦  <  A#1  > 

♦  ]  <  Ao  >  Gi  Hi  <FVlH‘  (27) 


where: 
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Hj  i  e  if ♦<*♦*>  *1°  (28) 

<  0,0  >  is  given  by  eq.  (25), 

<n6  the  remaining  symbols  are  os  defined  above. 
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For  smoothening  films,  our  result  is 
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where  <  ?  >  «  “f“  with  e  the  RMS  roughness  height  in 
each  inter  face,  ond  where  Y  ond  t  ore  os  previously 

defined. 

■We  note  that  the  effect  of  smoothening  films  is 
very  similar  to  that  of  interlayer  diffusion.  This  is 
particularly  true  if  the  RMS  roughness  height  is  at 
least  moderately  small,  so  that  *'•’  <  f  >  >  q  I  el, 
where  q  is  the  number  of  layers  within  one  longitudinal 
autocorrelation  length  of  the  roughness.  (The  analysis 
assumes  that  q|  *1  «  I  ). 

In  this  cose  the  intensity  of  the  diffuse  beam 
will  be  small  compared  to  that  of  the  specular  beam, 
even  though  the  specular  reflectivity  may  be 
considerably  less  than  it  would  be  in  the  absence  of 
roughness.  The  main  effect  of  the  roughness  in  this 
cose  is  to  decrease  the  reflectivity  of  the  multilayer, 
rather  than  to  scatter  radiation  diffusely.  This  is 
because  radiation  that  is  diffusely  scattered  from 
different  layers  odds  only  incoherently  under  the 
assumptions  of  the  smoothening  film  model. 

For  comparison,  suppose  one  models  the  diffusion 
process  by  considering  it  to  cause  the  ideal  multilayer 
profile  Atz)  to  be  convolved  with  some  diffusion  profile 
g(z). 

This  causes  the  parameter  y  to  be  multiplied  bv 
the  Fourier  transform  of  g(z).  If  g  is  taken  to  be  a 
Gaussian,  g(z)  =  exp(-.5(z/o)2),  the  reflected 
intensity  is  given  by  cr>  expression  identical  to  eq.(29) 
above  for  R,„„  in  the  presence  of  smoothening  films. 

Ref.  8  gives  an  expression  for  specular 
reflectivity  in  the  presence  of  identical  films  that 
applies  under  the  assumptions  of  the  scalar  model.  The 
same  egression  is  used  in  ref.  11. 

We  also  note  that  for  all  three  types  of 
roughness,  walkoff  considerations  can  be  shown  to  imply 
that  the  scalar  model  is  applicable  only  when  the 
seprration  between  the  specular  and  diffuse  beams  is 
within  the  acceptance  mgle  of  the  multilayer.  Since 
the  acceptance  angle  is  likely  to  be  of  the  order  of 
the  field  of  view  in  imaging  applications,  the  scalar 
model  can  be  opplied  to  m«iy  imaging  problems  of 
interest. 
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The  angle  of  maximum  reflectivity  from  a  multilayer  x-ray  reflector  is  influenced  by  absorption 
in  the  medium.  Using  multilayer  theory  we  show  how  the  full  refractive  correction  compares  to 
the  usual  correction  that  includes  only  the  effect  of  dispersion  The  physical  significance  of  the 
two  corrections  is  discussed.  The  absorption  correction  is  computed  for  some  multilayer  systems 
of  current  interest. 
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Much  progress  has  been  reported  recently  in  the  fabri¬ 
cation  of  artificial  Bragg  reflectors  for  the  soft  x-ray  re¬ 
gime  Applications  for  these  multilayer  systems  include 
monochromators  for  synchrotron  radiation,'  norma)  inci¬ 
dence  mirrors  for  subkilovolt  x-ray  imaging  systems."  graz¬ 
ing  incidence  mirrors  for  the  reflection  of  higher  energy  x 
rays  than  is  otherwise  practical,4  and  cavity  mirrors  for  fu¬ 
ture  x-ray  lasers 7  In  multilayers  used  with  such  devices,  the 
wavelength  [A ),  glancing  angle  of  incidence  10 ),  and  period 
spacing  (d  I  must  be  carefully  matched. 

It  is  well  know  that  the  simple  Bragg  relation 

2d  sin  0,  —  mA  !  1 ) 

must  be  corrected  for  refractive  effects.  The  crystal  plane 
spacing  d  and  wavelength  in  vacuum  A  are  both  measured  in 
angstroms.  0,  is  the  grazing  angle  of  incidence  and  the  inte¬ 
ger  m  is  the  order  of  diffraction.  The  standard  treatment  of 
the  refraction  correction  in  a  perfect  crystal  ignores  the  ef¬ 
fect  of  absorption.19  Expressed  in  angular  units,  this  stan¬ 
dard  correction  for  a  crystal  having  a  spatially  averaged  in¬ 
dex  of  refraction  n  =  1  —  $  —  i$  is 

0  —  0,  =  £/sin  dt  cos  dt ,  (2) 

where  it  is  assumed  that  6<  1  and  #  =  0. 

Equation  (2)  does  not  correctly  give  the  position  of 
maximum  reflectivity  in  a  system  when  n  is  complex.  Shack- 
lett  and  DuMond  have  suggested  that  the  absorption  correc¬ 
tion  be  made  graphically.10  Miller  has  derived  an  analytic 
expression  for  the  full  refractive  correction  using  the  Dar- 
win-Prins  dynamical  theory  of  crystal  diffraction."  In  this 
letter,  the  angular  position  of  diffraction  maxima  have  been 
obtained  by  simpler  means,  using  a  formalism  developed  to 
treat  multilayer  reflection. 

■JNt  note  that  in  many  applications  such  as  those  cited 
above,  the  incident  beam  can  be  treated  as  linearly  polarized. 
In  a  previous  report, 11  it  has  been  shown  that  in  this  case  the 
intensity  reflectance  R  of  a  multilayer  containing  a  large 
number  of  layers  is 


—  ID  +  iB) 


(3) 


t-Wt’-V  +  iB)2]  I  ' 

Here  p  is  the  amplitude  reflectivity  of  the  multilayer  and 
—  4 D  +  iff )  is  the  amplitude  reflectivity  of  a  single  priod  f  is 
a  detuning  parameter  defined  in  terms  of  the  optical  phase 


thickness  V  of  the  multilayer's  basic  period  through  the 
relation 


£  =  4*  —  mirm 
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where  we  have  defined 

^=—-1  dzn|z|(v/I  -  cos;  0/n:|z|] 

A  J  period 

m  2nd  sing  L  6  +  i0  \  (5| 

/l  \  sin:0  )' 

with  2  the  coordinate  perpendicular  to  the  layer  interfaces 
ThesecondpartofEq.  |5)follows because  n\z\  -  1 ,  <  1  in  the 
x-ray  regime. 

In  Ref  12,  the  parameters  D  ■+  iB  and  6  +  i$  are  de¬ 
rived  for  multilayers  consisting  of  alternating  layers  of  two 
materials  having  indicies  of  refraction  n ,  and  n3. 6  and  0  are^ 
the  averaged  unit  decrements  given  by  6  =  yS,  «+•  (1  —  yl  6: 
with  similar  expression  for#  where  y  is  the  division  param¬ 
eter  defined  by  ymd,/\dt  +  d3\.  D  +  iB  equals  -(-IT 
2r  sin  (mtry),  where  r  is  the  Fresnel  reflection  coefficient  To 
first  order  in  the  unit  decrement  differences,  r  is  given  by 

r  =  ((d<5  +  id#  1/2  sinJ  6  J P(0 ),  (6| 

where  P{8 1  =  1  or  cos  (20 1  depending  on  whether  the  inci¬ 
dent  beam  is  5  polarized  or  P  polarized,  respectively,  and 
d £  +  id#  *(5|  +  (#,)  —  (6-  ■+■  i0 j). 

We  now  calculate  the  positions  of  maximum  reflectiv¬ 
ity  D  and  B  are  virtually  independent  of  £.  since  6mBt.  so 
that  using  Eq.  (3), 

J_"  ^WRe(l^)«/Re,y-,. 

1R  de  \pd£de)  \  p  di  J 

<71 

where 

{2rsin(mfl7)]5j 

and 
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In  Eqs.  (8)  and  (3)  the  root  with  the  positive  imaginary  part 
must  be  chosen. 
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At  dB  /dB  m  0,  the  reflectivity  is  maximized.  There¬ 
fore,  the  petitions  of  the  diffraction  maxima  can  be  obtained 
from  the  condition 

ReF-0,  (9) 

where  we  have  used  the  identity  Re(F  " ')  =  |  Y  |-J  Re  (Y ). 

It  is  shown  in  Ref.  13  that  the  variable 

Q,~mir-Y  (10) 

is  the  quantity  known  as  the  equivalent  phase  thickness,  that 
has  been  introduced  into  multilayer  analysis  by  Herpin.14 

Equation  (9)  is,  therefore,  a  requirement  that  the  real 
part  of  the  multilayer's  equivalent  phase  thickness  equal  m  fl¬ 
at  the  position  of  maximum  reflectivity.  This  is  not  the  same 
as  the  conventional  requirement  that  the  real  pan  of  the 
optical  phase  thickness  'X  equal  mir. 

We  may  explain  the  difference  as  follows.  The  oscilla¬ 
tion  of  the  forward  traveling  wave  across  each  period  of  the 
multilayer  is  given  by  e  ~  **’,  since  <P,  is  the  equivalent  phase 
thickness.  This  oscillation  is  determined  not  only  by  the  opti¬ 
cal  phase  thickness  of  the  period  being  traversed,  but  also  by 
multiple  reflections  involving  the  other  periods  of  the 
multilayer. 

This  may  be  seen  by  eliminating  the  radical  in  the  de¬ 
nominator  of  Eq.  (3)  to  obtain 

-(£+  Y)/(D  +  iB),  (11) 

so  that  Y  *=  —  f -p  (D  +  iff )  and 

=  4>  +  p(D  +  iB),  (12) 

and  therefore 

,--'  =  e-1*[l-v>(D  +  /ff)]+0(r2).  (13) 

Here  the  term  e  ~  '*  represents  the  effect  of  transmission 
through  the  period.  The  term  +  iff  (represents  a 

multiple  reflection  process  consisting  of  transmission 
through  the  period,  reflection  from  all  succeeding  periods, 
and  reflection  back  into  the  initial  direction  by  the  period 
initially  traversed.  (Higher-order  multiple  reflections  are  of 
order  r2  and  are  small.)  A  phase  change  occurs  during  the 
multiple  reflection,  and  if  n  is  complex,  the  overall  oscilla¬ 
tion  of  the  field  across  the  period  is  shifted.  For  brevity's 
sake,  we  have  described  the  calculation  of  this  effect  as  a 
correction  for  absorption 

We  can  convert  the  resonance  condition  of  Eq.  (9)  to  a 
computationally  more  useful  form  as  follows:  The  condition 
Re  Y  *  0  implies  that  Y2  must  be  a  negative  real  number.  It 
is  readily  shown  that  the  condition  Im(  Y 2)  —  0  is  sufficient  to 
ensure  Ref  Y2)  < 0.  From  Eqs.  |S)  and  (6),  the  requirement 
that  ImfK2)  *  0  yields 


(sin  6  —  am  8,  )sia  9  ■  i  - 
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where  we  have  set  P {0  )»ff  (ff#).  Introducing  A0m8- 
and  using 

(sin  0- tin  0,)sin  0-40(004  0,  -  $40sin  0,)ain  0, 

+  Of402eoaJ0,) 
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(which  is  accurate  even  near  normal  incidence  where  the  two 
terms  in  the  first  parentheses  on  the  right  may  be  compara¬ 
ble),  we  find  that  40  is  given  by 

40  =  |1  -  [V\  -2»7sec20,)|/tan0,,  (16) 

where  q  is  defined  to  be  the  right-hand  side  of  Eq.  (14). 

Equation  (16)  is  valid  for  all  angles  outside  the  grazing 
incidence  regime,  except  that  no  physical  solution  exists  if 
the  radicand  on  the  right  side  is  negative,  which,  to  a  good 
approximation,  can  be  assumed  to  occur  if  0,  >  rtf 2 
-  \V2fi\- 

Away  from  normal  incidence  [i.e.,  it/ 2  —  0B>(V2 17)], 
Eq.  (16)  can  be  written  as 


$  A6A0sm2{mirr)P-{0g\ 

at)  =  -  —  - — — tt  - .  (17) 

sin  6t  cos  6t  m2irp  sin  6,  cos  6B 

The  first  term  is  the  usual  dispersion  correction.  The 
second  term  is  due  to  absorption. 

Equation  (16)  can  be  applied  to  a  general  noncentro- 
symmetric  structure  (such  as  a  crystal)  by  letting 
V  —  S  —  (DB  +  DB )  sin4  dt/(2m2tri0 ).  D,  ff,  D,  and  ff  are 
standard  parameter?  of  crystal  diffraction  theory ,  whose  use 
in  multilayer  analysis  is  discussed  in  Ref.  12. 

In  the  common  case  of  a  crystal  that  is  centrosymme- 
tric,  Eq.  (17)  may  be  applied  in  the  form 


sin  0,cos  6t 


[rJA /2ir]2F,F}P2(0g) 
B  sm0,  cos  6t 


(I81 


Here  F,  and  F2  are  the  real  and  imaginary  pans,  respectively, 
of  the  crystal  structure  factor.  r„  is  the  classical  electron  radi¬ 
us  and  d  is  the  reciprocal  of  the  unit  cell  volume. 

The  relative  magnitude  of  the  two  corrections  are  com¬ 
pared  in  Fig.  I  for  the  case  of  a  particular  reflecting  struc¬ 
ture.  The  ratio  plotted  is 


RATIO  =  (0,  -  02)/(0,  -  6B ),  (19) 

where  0:  is  the  true  Bragg  angle  as  corrected  for  absorption 
and  dispersion  and  0,  is  the  Bragg  angle  as  corrected  for 
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Refraction 

correction 

Absorption 

correction 

Ratio 

1*1 

Vanadium/ 
carbon 
(multilayer  *| 

1J910  3 

13.621  1 

0.263  8 

0  031  3 

-  10  61 

Tungsten/ 

carbon 

(multilayer*) 

14.317  4 

13.621  1 

0.869  1 

-0.1390 

1904 

TlAP 

(crystal) 

6S.826  37 

65.671  29 

0.164  28 

-  0.009  57 

6  19 

RAP  • 

(crystal! 

65.444  46 

65.283  96 

0.164  02 

-0003  80 

237 

KAP 

(crystal  l 

63.000  39 

62  859  98 

0.140  98 

-0.000  74 

5.28 

NAAP 

(crystal: 

63.421  85 

63.286  17 

0.135  73 

-0.000  19 

0  14 

nhap 

(crystal  i 

64  433  82 

64.304  00 

0  130  36 

-0.000  71 

0.55 

"Id  -  100  A  and  y  «  0  .5 


dispersion  only.  Equation  (19)  is  the  ratio  of  the  absorption 
correction  to  the  net  correction. 

To  obtain  Fig.  1,  we  have  used  preliminary  values  of 
atomic  scattering  factors  from  a  forthcoming  compilation. 15 

One  can  show  that  the  absorption  correction  must  al¬ 
ways  be  smaller  than  the  dispersion  correction,  but  the  plot 
shows  that  there  are  regimes  in  which  the  two  effects  are 
comparable  to  order  of  magnitude. 

First,  the  absorption  has  an  increased  effect  on  the  reso¬ 
nance  angle  near  normal  incidence,  where  the  correspon¬ 
dence  between  angular  shifts  and  shifts  in  phase  thickness  is 
nonlinear. 

Also,  the  absorption  correction  is  comparatively  large 
at  wavelengths  just  above  the  carbon  K  edge,  and  a  short 
wavelengths;  in  these  regions  the  effect  of  multiple  reflec¬ 
tions  is  comparatively  strong. 

Table  1  shows  the  refraction  and  absorption  corrections 
associated  with  reflection  from  other  periodic  structures.  In 
all  cases  m  *  1  and  A  m  23.6  A  (oxygen  Ka  line).  62  is  the 
predicted  diffraction  peak  position,  valid  to  all  orders  of  A  6, 
computed  using  Eq  ( 14).  The  refraction  and  absorption  cor¬ 
rections  are  the  first  and  second  terms  of  the  right-hand  side 
of  Eq.  (17),  respectively.  The  difference  between  the  82  col¬ 
umn  and  the  sum  of  the  next  three  columns  reflects  the  con¬ 
tribution  of  higher  order  A$  terms  which  have  been  ignored 
in  Eq.  (17). 

Finally,  the  last  column  is  the  ratio,  Eq.  (19),  in  percent 
of  absorption  correction  to  the  net  correction.  As  long  as  the 
refractive  indices  have  real  pails  less  than  I,  the  refraction 
correction  always  adds  to  the  Bragg  angle  9, .  The  absorp¬ 
tion  correction  generally  subtracts  from  0*;  however,  there 
are  cam  where  absorption  adds  to  9,  as  is  evident  in  the  case 
of  the  vanadium/carbon  multilayer. 
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